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Mid-to-High-Frequency  Bottom  Loss 
in  the  East  China  Sea 

Jee  Woong  Choi  and  Peter  H.  Dahl 


Abstract — Bottom-loss  measurements  made  in  the  East  China 
Sea  in  May-June  2001  as  part  of  the  Asian  Sea  International 
Acoustics  Experiment  as  a  function  of  frequency  (2-20  kHz)  and 
seabed  grazing  angle  (15°-24°)  are  presented.  The  measurements 
are  interpreted  as  estimates  of  the  modulus  of  the  plane  wave 
reflection  coefficient  and  data  are  compared  to  predicted  values 
using  a  reflection  coefficient  model,  based  on  a  two-layered  sedi¬ 
ment  for  which  the  sound  speed  in  the  surficial  sediment  layer  is 
allowed  to  vary  as  a  linear  k 2  profile,  where  k  is  acoustic  wave 
number.  The  region  below  this  layer  is  modeled  as  a  half-space 
with  constant  density  and  sound  speed.  The  reflection  coefficient 
model  is  driven  by  eight  geoacoustic  parameters;  these  are  esti¬ 
mated  from  the  data  by  minimizing  the  weighted  squared  error 
between  the  data  and  the  model  predictions  for  a  candidate  set  of 
parameters.  The  parameter  estimates  for  the  sediment  layer  are 
thickness,  0.9  ±  0.5  m;  density,  2.0  ±  0.1  g/cm3;  and  attenu¬ 
ation,  0.25  ±  0.05  dB/m/kHz,  with  sediment  layer  sound  speed 
increasing  from  1557  ±  4  m/s  at  the  water-sediment  interface 
to  1625  ±  35  m/s  at  a  depth  of  0.9  m.  The  parameter  estimates 
for  the  half-space  are  density,  2.0  ±  0.1  g/cm3;  attenuation, 
0.25  ±  0.15  dB/m/kHz;  and  sound  speed,  1635  ±  52  m/s. 
Variances  for  these  estimates  are  derived  using  the  Bootstrap 
method.  This  parameter  set  produced  model  curves  that  agreed 
reasonably  well  with  the  observations  of  bottom  loss  over  the  entire 
frequency  range  and  is  consistent  with  the  range  of  independently 
measured  geoacoustic  variables.  Since  this  mid-to-high-frequency 
data  set  does  not  provide  detailed  information  about  the  sediment 
structure  for  depths  beyond  about  3  m,  the  geoacoustic  parameter 
set  is  more  properly  viewed  as  description  of  the  sediment  layer 
and  sediments  in  the  underlying  2  m.  Similarly,  a  self-consistent 
construction  of  a  geoacoustic  model  for  the  East  China  Sea  should 
necessarily  amalgamate  the  mid-to-high-frequency  results  given 
here  with  results  obtained  at  lower  frequencies. 

Index  Terms — Bottom  loss,  bottom  reflection,  geoacoustic  pa¬ 
rameters,  ocean  acoustic  propagation,  sediment  layering. 


I.  Introduction 

THIS  paper  presents  results  of  an  experiment  to  measure  the 
modulus  of  the  mid-to-high-frequency  bottom  reflection 
coefficient,  via  bottom-loss  measurements,  at  one  site  in  the  East 
China  Sea  as  part  of  the  2001  Asian  Seas  International  Acoustics 
Experiment  (ASIAEX).  The  measurements  were  conducted  350 
km  east  of  Shanghai,  China,  in  the  Chinese  continental  margin 
(Fig.  1),  at  location  29°  39’N,  126°  49’E  in  waters  105  m  in 
depth,  in  an  area  historically  characterized  by  very  fine  sand 
sediment  [1]. 

Manuscript  received  June  11,  2003;  revised  September  3,  2003.  This  work 
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Fig.  1.  Location  of  the  ASIAEX  East  China  Sea  experimental  site.  The 
measurements  were  made  within  O  (100)  m  of  the  position  29°  39’N,  126° 
49’E;  water  depth  is  105  m. 

Data  were  collected  as  function  of  frequency  (range  2-20 
kHz)  and  grazing  angle  (range  15°-24°)  and  the  observable  is 
an  estimate  of  bottom  loss  based  on  averaged  intensity.  This  ob¬ 
servable  is  compared  to  -20  log  \R\,  where  \R\  is  the  modulus 
of  the  plane  wave  reflection  coefficient.  The  range  of  grazing 
angle,  though  limited,  was  designed  to  span  the  expected  crit¬ 
ical  angle,  which  is  a  key  factor  influencing  sound  propagation 
in  a  shallow-water  waveguide  and  a  point  about  which  \R\  will 
show  stronger  variation  with  grazing  angle.  Observations  are 
compared  to  predictions  for  which  the  forward  problem  con¬ 
sists  of  computing  \R\  associated  with  a  two-layered  sediment 
underlying  a  water  layer,  based  on  candidate  sets  of  geoacoustic 
parameters.  An  optimum  set  of  geoacoustic  model  parameters  is 
identified  from  an  exhaustive  search  over  a  constrained  param¬ 
eter  space  based  on  a  weighted  least-squared  error  criterion. 

Geoacoustic  ground  truth  measurements  of  the  experimental 
site  were  made  in  2000  and  2001  as  part  of  the  ASIAEX  East 
China  Sea  field  program.  These  are  in  the  form  of  gravity  and 
piston  cores  and  water-gun  and  chirp  sonar  generated  subbottom 
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Antenna,  battery  pack 


Fig.  2.  Acoustic  source  system  (left)  and  receiving  vertical  line  array  system  (right)  used  to  measure  bottom  loss.  The  data  presented  in  this  paper  are  from  a 
single  hydrophone  from  each  of  the  two  vertical  line  arrays. 


profiles  (see  Miller  et  al.  [1]  for  detailed  discussion  of  geoa¬ 
coustic  measurements).  In  short,  these  data  suggest  a  subbottom 
structure  at  the  ASIAEX  site  that  consists  of  a  layer  of  sedi¬ 
ment  of  variable  thickness  starting  at  the  water-sediment  inter¬ 
face,  for  which  near  the  site  of  acoustic  measurements  the  thick¬ 
ness  was  of  O  (1  m).  Coring  data  [1]  obtained  within  100  m  of 
the  site  of  acoustic  measurements  indicate  a  sound  speed  going 
from  approximately  1575  m/s  at  the  water-sediment  interface  to 
1 600- 1 675  m/s  at  approximately  1  m  into  the  sediment.  There  is 
also  evidence  of  centimeter-scale  layering  superimposed  on  this 
gradual  increase.  The  coring  data  also  show  a  mean  grain  size 
(in  phi  scale)  of  3.470,  consistent  with  very  fine  sand,  and  bulk 
density  between  1.85-2.05  g/cm3.  The  root  mean  square  (rms) 
bottom  relief  was  0.005  m  within  a  4-m  measurement  aperture 
based  on  the  measurements  of  fine-scale  bottom  properties  by 
Tang  [2],  made  at  the  precise  location  of  the  acoustic  bottom 
loss  measurements  reported  here. 

This  paper  is  organized  as  follows.  The  field-measurement 
procedure  and  estimation  of  bottom  loss  are  described  in  Sec¬ 
tion  II  and  results  are  presented  in  Section  III.  Section  IV  out¬ 
lines  the  bottom-reflection  model  based  on  a  two-layered  fluid 
sediment  for  which  a  sediment  layer  with  positive  sound  speed 
gradient  is  sandwiched  between  sea  water  and  a  homogeneous 
half-space.  This  model  is  used  in  forward  calculations  of  the  re¬ 
flection  coefficient,  which  in  turn  are  used  in  Section  V  in  an  ex¬ 
haustive  search  for  geoacoustic  parameters  that  minimize  error 
between  modeled  and  measured  bottom  loss.  A  summary  and 
the  conclusion  are  given  in  Section  VI. 

II.  Field  Measurments  and  Estimation  of  Bottom  Loss 

Bottom-reflection  loss  measurements  were  made  on  May  3 1 
and  June  5, 2001  and  were  part  of  set  of  nominally  0.5-km  range 
propagation  measurements  designed  to  study  both  sea  surface 


and  seabed  scattering  and  reflection.  Some  results  pertaining  to 
the  sea  surface  are  discussed  in  [3].  An  acoustic  source  system 
was  deployed  off  the  stem  of  the  R/V  Melville  (while  it  main¬ 
tained  dynamic  positioning)  at  a  depth  of  either  25  or  50  m. 
A  family  of  five  transmit  waveforms,  CW  pulses  with  a  length 
of  2  or  3  ms  and  center  frequency  equal  to  2,  4,  8,  16,  and  20 
kHz,  plus  two  additional  waveforms,  was  cycled  every  100  s  for 
20  times;  data  corresponding  to  a  set  of  20  x  7  transmissions 
was  identified  as  an  experimental  set,  with  this  process  repeated 
after  a  few  seconds  delay.  An  ITC  2010  transducer  from  Inter¬ 
national  Transducer  Corporation,  Santa  Barbara,  CA,  was  used 
for  the  2-kHz  transmission  and  an  ITC  (omnidirectional)  1007 
transducer  was  used  for  all  other  transmissions.  At  2  kHz,  the 
transmit  pattern  for  the  ITC  2010  used  in  this  study  is  invari¬ 
able  within  the  vertical  angular  range  0°-20°  and  there  is  ap¬ 
proximately  a  0.5-dB  change  between  20°-25°  at  most.  Thus, 
both  transmitting  systems  are  considered  to  be  omnidirectional 
over  the  angular  range  used  in  the  measurements.  The  trans¬ 
missions  were  received  on  two  colocated  vertical  line  arrays, 
each  of  length  1  m,  with  one  array  at  depth  26  m  and  the  other 
array  at  depth  52  m,  referenced  to  the  array  center  (Fig.  2).  The 
two-array  system,  called  moored  receiving  array  (MORAY), 
was  deployed  at  a  nominal  range  of  500  m  from  the  source  with 
precise  range  depending  on  measurement  time  and  established 
acoustically.  Data  for  this  analysis  originate  from  single  omni¬ 
directional  elements  (ITC  1042)  on  the  26-  and  52-m-depth  re¬ 
ceiving  array.  The  data  were  sampled  at  50  kHz  and  sent  back 
to  the  Melville  through  a  radio-frequency  (RF)  modem. 

Sound  speed  profiles  were  monitored  with  frequent  con¬ 
ductivity-temperature-depth  (CTD)  casts  made  from  both  the 
Melville  and  the  nearby  Chinese  R/V  Shi  Yan  3.  In  general,  the 
sound  speed  versus  depth  variation  in  the  East  China  Sea  is 
linked  to  the  tides  during  the  late  spring  and  summer  months, 
and  the  sound  speed  in  the  upper  mixed  layer  (nominal  depth 
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Fig.  3.  (a)  First  five  eigenrays  based  on  the  sound  speed  profile  for  set 

10.  Source  depth  and  range  between  source  and  receiver  are  25  and  499  m, 
respectively,  (b)  Receiver  depths  26  m  for  the  upper  array  and  (c)  52  m  for  the 
lower  array. 


30  m)  often  varied  between  1525-1530  m/s,  while  the  sound 
speed  near  the  water-sediment  interface  was  a  stable  1519  m/s. 

Fig.  3  shows  propagation  conditions  for  measurement  set  10, 
obtained  at  1030  UTC  May  31,  and  for  which  the  source  depth 
was  set  to  25  m;  (a)  shows  the  sound  speed  profile;  (b)  and  (c) 
show  eigenrays  corresponding  to  the  first  five  arrivals  (direct, 
surface,  bottom,  surface-bottom,  and  bottom-surface),  based  on 
the  sound  speed  profile  in  (a)  and  computed  using  a  generic  ray- 
based  propagation  code.  The  upper  and  lower  receiving  arrays 
each  sample  slightly  different  bottom  bounce  grazing  angles, 
in  this  case  18.1°  for  the  upper  array  and  15.3°  for  the  lower 
array.  A  range  of  bottom-grazing  angles  between  15°-21°  was 
achieved  through  using  the  combination  of  upper  and  lower  re¬ 
ceiving  arrays,  the  two  source  depths  at  25  and  50  m,  plus  varia¬ 
tion  in  source  to  receiver  range  over  the  course  of  the  experiment 
that  was  between  400-500  m.  Finally,  this  grazing  angle  range 
was  extended  to  24°  by  exploiting  data  that  involved  a  single 
sea-surface  interacting  path  in  addition  to  a  bottom  bounce  path. 
Clearly,  coherent  reflection  from  the  sea  surface,  for  which  the 
rms  wave  height  varied  but  was  never  below  about  0.2  m,  was 
negligible.  Instead,  intensity  is  all  incoherent;  sea-surface  in¬ 
teraction  imparts  no  energy  loss  provided  all  intensity,  which 
undergoes  both  time  and  angle  spreading,  is  measured.  This  cri¬ 
terion  is  satisfied  with  our  wide-beam  widths  and  receive  time 
window.  Energy  loss  associated  with  attenuation  from  near  sur¬ 
face  bubbles  can  occur,  but  this  effect  is  negligible  under  exper¬ 
imental  field  conditions  sampled  for  these  measurements  char¬ 
acterized  by  mild  wind  speed  (less  than  about  5  m/s)  [4],  [5]. 
Thus,  in  these  cases  the  measured  loss  in  bottom-interacting 
paths  in  excess  of  that  due  to  chemical  absorption  and  geo¬ 
metric  spreading  (both  of  which  must  be  modeled),  is  taken  to 
be  associated  with  the  bottom  loss.  The  bottom- loss  estimates  at 


Time  relative  to  Direct  Path  (ms) 


Fig.  4.  Individual  pings  (thin  lines)  and  average  intensity  based  on  a  20-ping 
average  (thick  line)  of  16-kHz  data  for  the  lower  receiver  array  for  set  10. 


grazing  angles  greater  than  21°,  however,  are  encumbered  with 
a  higher  variance  owing  to  fluctuations  imparted  by  interaction 
with  the  sea  surface. 

Fig.  4  shows  the  received  intensity  level  for  20  individual 
pings  (thin  lines)  and  their  average  intensity  (thick  line)  for  the 
16-kHz  measurements  received  on  the  lower  receiver  array  for 
set  10.  The  intensity  level  is  referenced  to  that  of  the  mean  peak 
level  for  the  direct  path,  which  here  is  assigned  (arbitrarily  in 
this  figure)  to  be  0  dB.  The  ray-based  spreading  and  absorp¬ 
tion  loss  estimate  for  the  direct  path  is  57.1  dB  and  that  for  the 
bottom  bounce  path  is  54.7  dB.  In  this  case,  the  direct  path  has 
slightly  higher  spreading  loss  owing  to  refraction  in  the  water 
column.  Absorption  loss  was  predicted  with  [6]  using  salinity 
and  temperature  measured  by  CTD.  These  values,  in  combina¬ 
tion  with  the  mean  level  for  the  bottom-interacting  path,  -2.5 
dB,  give  an  estimate  for  -20  log  |/?|  equal  to  4.9  dB  at  15°  for 
16  kHz.  However,  a  total  energy  measure,  based  on  an  integral 
over  a  time  window  somewhat  greater  than  the  pulse  length,  is 
used  to  bottom  loss  in  the  manner  of  [7].  This  estimate  puts 
—20  log10  |/?|  equal  to  4.4  dB,  which  becomes  our  reported  esti¬ 
mate.  Statistical  fluctuations  in  the  level  of  the  single  bottom-re¬ 
flected  path  based  on  20-ping  average  are  ±1.4  dB,  which  are 
added  in  quadrature  to  the  estimated  error  for  source  level,  ±0.9 
dB,  giving  a  total  uncertainty  for  this  estimate  as  ±  1 .7  dB.  Note 
that  total  measurement  uncertainty  varied  depending  on  mea¬ 
surement  time  of  day  and  frequency.  For  grazing  angles  for  es¬ 
timates  involving  the  sea-surface  bounce  path,  total  uncertainty 
was  dominated  by  statistical  fluctuations,  which  put  this  value 
closer  to  ±3  dB.  Finally,  the  total  energy  measure  is  essential 
for  data  that  involved  a  single  sea-surface  interacting  path  be¬ 
cause  of  time  spread.  In  this  case,  the  length  of  time  window 
must  be  wide  enough  to  contain  the  effect  of  time  spreading  of 
the  acoustic  signal  caused  by  sea  surface  roughness. 

For  the  single  bottom-reflection  path,  minimal  pulse  elonga¬ 
tion  is  observed,  owing  that  we  believe  to  the  bottom  boundary 
being  acoustically  smooth,  as  evidenced  by  Tang’s  [2]  estimate 
for  rms  seabed  relief  H  being  0.005  m.  A  measure  of  rough¬ 
ness  is  the  Rayleigh  parameter  \  =  2  kH  sin(09),  where  k 
is  acoustic  wave  number  and  6g  is  the  nominal  grazing  angle. 
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Fig.  5.  Measured  bottom  loss  as  a  function  of  grazing  angle  and  frequency 
compared  to  model  for  —  201og|/?|  based  on  a  homogenous  half-space 
representation  for  the  seabed. 

At  16  kHz  and  grazing  angle  15.3°  \  =  0-18  and  in  general 
X  <  0.34  for  all  measurements  based  on  the  highest  frequency 
and  grazing  angle,  20  kHz  and  24°,  respectively,  indicating  that 
the  bottom  is  an  acoustically  flat  surface  and  the  reflected  field 
is  primarily  generated  within  the  first  Fresnel  zone  [8],  [9].  At  2 
kHz,  the  first  Fresnel  zone  for  the  single  bottom-reflected  signal 
for  the  upper  receiver  for  set  10  is  an  approximate  ellipse  61 
m  x  19  m  and  at  20  kHz  this  ellipse  is  reduced  to  19  m  x  6 
m;  for  the  lower  receiver,  the  Fresnel  zone  is  slightly  larger. 
The  Fresnel  zone  dimensions  for  the  surface-bottom  path  are 
of  the  same  order.  Thus,  these  measurements  represent  an  in¬ 
herent  range  average  over  ranges  of  order  10-100  m. 

III.  Measurement  Results 

Fig.  5  shows  estimates  of  bottom  loss  expressed  in  decibels 
as  a  function  of  grazing  angle  and  frequency  and  assumed  to 
correspond  to  -20  log  \R\.  Measurements  for  which  the  grazing 
angle  is  greater  than  21°  are  based  on  paths  that  interacted 
once  with  the  sea  surface  and  have  greater  uncertainty  owing 
to  fluctuations  imparted  by  interaction  with  the  sea  surface. 
Measurement  variance  at  all  grazing  angles  tends  to  increase 
with  frequency.  The  scintillation  index  (normalized  variance 


of  intensity)  for  the  bottom  reflected  path  generally  goes  from 
about  0.01  to  0.2  for  the  frequency  range  2-20  kHz,  increasing 
roughly  linearly  with  frequency.  Therefore,  we  assume  that 
the  increasing  spread  of  measured  data  with  frequency  is 
associated  with  fluctuations  linked  to  small-scale  sound  speed 
changes  in  the  water  column.  Other  sources  of  variation  in 
the  data  are  postulated  to  be  associated  with  the  existence  of 
finer  (centimeter-scale)  layering  and  scattering  from  volume 
inhomogeneities  within  the  sediment;  however,  further  study  of 
these  effects  is  beyond  the  scope  of  this  paper. 

It  is  of  interest,  and  useful  for  reference,  to  first  compare 
our  measurements  to  a  homogeneous  fluid  half-space  model  for 
-20  log  \R\,  as  shown  by  the  curves  in  Fig.  5.  The  density,  sound 
speed,  and  attenuation  coefficient  for  this  model  are  1 .78  g/cm3, 
1600  m/s,  and  0.15  dB/m/kHz,  respectively.  These  geoacoustic 
parameters  were  estimated  by  Tang  [10]  using  ambient  noise 
within  the  1-3-kHz  band  recorded  on  a  vertical  line  array  that 
was  also  deployed  from  the  R/V  Melville.  The  half-space  model 
clearly  does  not  predict  the  variation  as  a  function  of  grazing 
angle  seen  in  the  measurements  at  frequencies  greater  than  2 
kHz.  Still,  the  half-space  model  does  reasonably  well  in  pre¬ 
dicting  the  2-kHz  data  (with  the  exception  of  the  two  measure¬ 
ments  made  at  grazing  angles  greater  than  22°). 

In  the  next  section,  sediment  layering  will  be  considered  as 
the  cause  of  frequency  and  angle  dependence  in  estimates  of 
bottom  loss,  and  a  layer  will  be  incorporated  into  a  forward 
model  for  the  seabed  reflection  coefficient.  Holland  [11]  has  also 
reported  similar  angular  dependence  in  bottom-loss  measure¬ 
ments  over  this  frequency  range,  his  measurements  illustrating 
an  influence  of  fine-scale  layering  on  reflection. 


IV.  Forward  Model  for  the 
Bottom-Reflection  Coefficient 

As  noted  previously,  the  subbottom  structure  in  the  East 
China  Sea  is  considered,  according  to  the  findings  of  Miller  et 
al.  [1],  to  have  an  O  (1-m)  thick  layer  of  sediment  starting  at 
the  water-sediment  interface  at  the  experimental  site.  Given 
this  finding,  we  postulate  a  forward  model  for  the  reflection 
coefficient  to  consist  of  a  water  layer  overlying  a  two-layered 
fluid  sediment.  This  can  be  expressed  as  [8] 


R  _  ^12(1  +  R23)  4-  R‘23 ( 1  +  R21) 
1  -  R21R23 


(1) 


where  Rij  is  the  reflection  coefficient  at  boundary  ij  and  sub¬ 
scripts  1-3  indicate  water,  finite-depth  layer,  and  half-space,  re¬ 
spectively.  Additional  specifications  of  the  forward  model  are 
as  follows.  Layer  2  (or  the  sediment  layer),  of  thickness  L,  is  al¬ 
lowed  to  have  some  positive  sound  speed  gradient,  consistent 
with  the  existence  of  sandy  sediments  [12].  Layer  3  is  mod¬ 
eled  as  an  infinite  half-space  of  constant  sound  speed  that  ter¬ 
minates  the  bottom  impedance.  Given  the  2-20-kHz  frequency 
range  and  <25°  grazing  angle  range,  moderate  sediment  atten¬ 
uation  renders  our  measurements  insensitive  to  sediment  prop¬ 
erties  greater  than  a  few  meters  in  depth;  therefore,  as  a  prac¬ 
tical  matter,  our  inversion  goal  is  to  identify  a  starting  sound 
speed  for  layer  3.  In  other  words,  we  have  no  information  in  our 
acoustic  data  to  justify  additional  details  in  layer  3,  e.g.,  such 
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as  additional  layered  structure  and  sound  speed  gradients,  etc. 
The  a  priori  information  that  the  sediment  thickness  L  is  O  (1 
m),  puts  kL  >  8  for  the  lowest  frequency  of  2  kHz.  Robins 
[13]  suggests  that  density  gradients  have  decreasing  influence 
on  bottom  loss  for  this  condition.  Thus,  the  density  in  layer  2 
is  assumed  to  be  some  constant  value,  with  same  assumption 
for  the  density  in  layer  3  (but  it  is  allowed  to  differ  from  that  in 
layer  2).  Finally,  constant  values  (but  not  necessarily  the  same) 
for  sediment  attenuation  are  assigned  to  layers  2  and  3,  and  sed¬ 
iment  attenuation  is  included  by  making  the  sound  speed  com¬ 
plex  in  the  manner  of  [  14].  We  believe  that  such  a  model  for  the 
bottom-reflection  coefficient  roughly  matches  the  inherent  re¬ 
solving  capability  of  the  measurements.  For  example,  looking 
ahead  to  the  next  section,  measurements  made  at  4  and  8  kHz 
clearly  exhibit  oscillations  as  a  function  of  grazing  angle  that 
can  only  be  reconciled  with  a  model  that  includes  some  kind  of 
layered  structure  in  the  upper  sediments. 

The  model  of  Ainslie  [15]  (see  also  [16])  accommodates 
the  previous  formulation  for  the  plane  wave  reflection  coeffi¬ 
cient  by  providing  analytical  solutions  of  the  depth-dependent 
Helmholtz  equation  for  certain  canonical  gradients.  Elements 
of  this  model  are  laid  out  as  follows.  The  boundary  condition 
that  the  pressure  and  normal  component  of  particle  velocity 
are  continuous  across  a  boundary  is  imposed  to  obtain  the 
reflection  coefficient  for  each  layer.  At  the  water-sediment 
layer  interface,  the  boundary  conditions  yield 


1  +  R\2 

—  (I-R12) 
Pi 


=  TnP+(z  =  0) 

=  riaa/»+(o) 

Pi  dz 


(2) 

(3) 


where  TtJ  is  the  transmission  coefficient  at  boundary  ij  and 
pt  and  7 i  are  sediment  density  and  vertical  wave  number  in 
layer  z,  respectively,  P2(z)  and  P2(z)  represent  the  down-  and 
up-traveling  fields  in  the  sediment  layer  and  the  substitution  of 
(2)  into  (3)  gives 


The  substitution  of  (7)  into  (8)  gives 


/?21 


/~(0)-Cl2 
/+( 0)  +  C12 


Using  the  similar  method,  R2 3  is  given  by 


Pfw  /+(/*) -C32 

P2-(/l)/-(/l)+C32* 


(9) 


(10) 


If  the  pressure  field  is  constant  with  depth,  f±(z)  =  1  and 
(6),  (9),  and  (10)  become  reflection  coefficients  for  layers  of 
homogeneous  fluid  media. 

Finally,  our  forward  model  with  which  to  match  the  observed 
data  (bottom  loss)  is  a  model  for  the  bottom-reflection  coeffi¬ 
cient,  for  which  the  sound  speed  profile  in  the  sediment  layer 
is  allowed  to  assume  a  linear  k2  profile.  For  this,  the  exact  ana¬ 
lytical  and  Wenzel,  Kramers,  and  Brillouin  (WKB)  approxima¬ 
tion  solutions  from  [15]  are  utilized.  This  forward  model  will 
be  governed  by  eight  parameters,  which  are  discussed  subse¬ 
quently.  Note,  however,  that  for  the  linear  k2  sound  speed  pro¬ 
file,  a  numerical  error  develops  when  the  sound  speed  gradient  is 
small  compared  to  acoustic  frequency.  This  error  can  be  signifi¬ 
cant  for  the  combination  of  high  frequencies  and  parameters  that 
describe  low  sound  speed  gradients.  To  include  such  gradients 
in  the  parameter  space  from  which  an  optimal  solution  will  be 
sought,  the  WKB  approximation  solution  is  used  in  the  forward 
model  for  the  field  in  the  sediment  layer  instead  of  the  exact  so¬ 
lution  when  deemed  necessary;  otherwise,  the  exact  solution  is 
employed.  The  known  WKB  failure  at  turning  points  is  not  en¬ 
countered  with  this  protocol  because  turning  points  within  the 
sediment  layer  do  not  arise  for  profiles  for  which  this  approxi¬ 
mation  is  applied. 


V.  Inversion  Based  on  a  Two-Layered 
Sediment-Reflection  Model 


P27i  dPf(o)/dz 

ppn  il2p+(o) 

mi  dPljoydz  • 

P 1T2  z72P+(0) 


(4) 


Ainslie  [15]  defines  the  right  terms  of  the  denominator  and  nu¬ 
merator  of (4)  as 


/*(*) 


dFfjdz 


(5) 


Defining  as  the  impedance  ratios  {pjli/ Pilj),  (4)  then  be¬ 
comes 


C12  -  /+ 

C12  +  /+(0) 


(6) 


Considering  the  reflection  (R2 1 )  of  the  up-traveling  signal  at  the 
water-sediment  layer  interface,  the  boundary  conditions  yield 


P2-(0)  +  /?2iP2+(0)  =  T21  (7) 

*'7i  Tn  1  pPf(O)  „  dP+(0)l 

- —  — ~ -  +  rC2i -  .  (8) 

Pl  p2  I  OZ  OZ 


Geoacoustic  parameters  for  the  reflection  coefficient  model  in 
Section  IV  are  determined  by  a  weighted  least-squared  error  cri¬ 
terion  for  the  total  mean-squared  error  between  a  model  vector 
based  on  candidate  set  of  parameters  and  the  data  vector.  Fixed 
values  are  used  for  the  sea-water  sound  speed  just  above  the 
sediment  interface  (1519  m/s)  and  density  ( 1 .0  g/cm3 )  based  on 
CTD  measurements,  while  an  eight-parameter  set  for  computing 
the  reflection  coefficient  is  varied  within  ranges  constrained  by 
a  priori  information  (Table  I).  This  parameter  set  consists  of:  1 ) 
layer  thickness  L\  2)  sound  speed  at  top  of  the  sediment  layer 
c2(0);  3)  density  of  the  sediment  layer  p2\  4)  attenuation  in 
the  sediment  layer  o2//  expressed  in  decibels/meter/kilohertz; 
5)  sound  speed  at  the  bottom  of  the  sediment  layer  c2(L)\  6) 
sound  speed  in  the  half-space  c3;  7)  density  of  the  half-space 
p3;  and  8)  attenuation  in  the  half-space  a3 //.  The  quantity 
Q  =  {Xi  -  Mt)TWe(Xt  —  Mi)  is  the  objective  function  where, 
X{  is  the  measured  bottom  loss  expressed  in  decibels,  M,  is 
model  prediction,  and  We  is  a  diagonal  weight  matrix  com¬ 
posed  of  l/a2  where  a2  is  the  measurement  variance  for  the 
nth  bottom  loss  estimate  (total  =175).  As  noted  previously,  data 
variance  depends  somewhat  on  grazing  angle  and  frequency  and 
the  objective  function  takes  into  account  the  data  uncertainty;  an 
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TABLE  I 

Constrained  Ranges  of  Geoacoustic 
Parameters  Used  for  the  Inversion 


Parameter 

Initial  search  space 

Final  search  space 

cj(  0) 

1520  1700  (m/s) 

1551  -  1562  (m/s) 

Ci(L) 

c?(0)  -  2000  (m/s) 

1600  -  1695  (m/s) 

c i 

1550-  2000  (m/s) 

1580-  1690  (m/s) 

P2 

1.0  -2.5  (g/cm3) 

1 .6  -  2.3  (g/cm') 

1.0  -  2.5  (g/cm') 

1.6 -2.3  (g/cm') 

L 

0.1  -  3.0  (m) 

0.4  -  1.5  (m) 

a  2/f 

0.1  -  0.7  (dB/m/kHz) 

0.1  -0.4  (dB/m/kHz) 

Otj/f 

0. 1  -  0.7  (dB/m/kHz.) 

0. 1  -  0.4  (dB/m/kHz) 

Fig.  6.  Diagram  of  the  best-fit  geoacoustic  parameters  for  two-layered 
sediment  model  and  associated  uncertainties. 

uncorrelated  error  assumption  is  invoked  insofar  as  constructing 
the  We  matrix. 

There  are  two  ranges  for  the  parameter  set  shown  in  Table  I. 
The  larger  initial  search  space  is  first  examined  using  a  coarser 
grid  with  which  to  discretize  each  parameter  range.  The  results 
of  this  first  stage  of  analysis  readily  identify  a  more  narrowed 
range  for  each  parameter,  which  becomes  the  final  search  space. 
It  is  this  space  over  which  an  exhaustive  search  is  completed 
to  minimize  the  objective  function  and  from  which  the  best-fit 
geoacoustic  parameter  set  is  selected. 

Results  of  the  inversion  for  the  best-fit  geoacoustic  param¬ 
eter  set  describing  a  sediment  layer  and  underlying  sediment 
half-space  are  diagrammed  in  Fig.  6.  Note  that  the  best  esti¬ 
mates  for  density  (2.0  g/cm3)  and  attenuation  (0.25  dB/m/kHz) 
in  the  sediment  layer  and  half-space  are  the  same,  which  re¬ 
flects,  in  part,  an  inability  to  resolve  additional  details  charac¬ 
terizing  geoacoustic  parameters  below  the  sediment  layer  with 
this  mid-to-high-frequency  data  set.  We  estimate  that  the  inver¬ 
sion  results  apply  to  depths  no  greater  than  about  2  m  into  the 


half-space  or  about  3  m  into  the  sediment  overall,  given  that  the 
lowest  frequency  is  2  kHz,  and  0.25  dB/m/kHz  holds  as  a  rea¬ 
sonable  estimate  for  attenuation.  The  data  provide  little  infor¬ 
mation  on  the  sediments  beyond  this  depth,  which  is  the  reason 
for  the  dashed  line  in  Fig.  6. 

An  additional  comment  on  attenuation  is  as  follows:  the 
relation  between  the  attenuation  in  decibels/meter  and  fre¬ 
quency  (/)  can  be  expressed  as  a  =  6/m,  where  b  is  a 
constant  and  estimates  of  m  have  been  reported  within  the 
range  of  0.5-2  [20].  The  attenuation  estimated  in  this  work, 
0.25  decibels/meter/kilohertz,  is  by  definition  equivalent  to  a 
linear  relation  (m  =1)  with  frequency.  This  result  should  not 
be  viewed  as  a  contradiction  of  any  proposed  nonlinear  trends, 
as  we  believe  that  the  one-decade  frequency  range  (2-20  kHz) 
in  our  data  set  is  not  large  enough  to  ascertain  a  frequency 
dependence  strongly  different  from  linear.  The  attenuation  is 
considered  by  letting  the  acoustic  wave  number  in  the  seabed 
be  complex.  The  complex  wave  number  ( k  =  k(l  +  id))  is 
specified  by  the  loss  tangent  6,  which  is  related  to  the  attenu¬ 
ation  coefficient  according  to 


which  implies  that  the  loss  tangent  is  depth  dependent,  but  the 
value  of  b  =  a/ f  is  a  constant  with  depth. 

The  linear  k2  sound  speed  profile  governing  the  reflection 
coefficient  model  is  given  by 

c2(z)  =  C2(0)[1  -  20z/c2(O)]~1/2  (12) 

where  ft  is  the  gradient  of  sound  speed,  which  can  be  approxi¬ 
mated  by 

c2(z)  «  c2(0)[l  +  [iz/c2(0)\  (13) 

for  |/J|z/c2(0)  <  1.  In  this  result,  fi  «  74.77  -  0.53i  (s-1), 
which  puts  \0\z/c2(O)  ~  0.04  at  z  =  L\  the  reason  for  the 
quasilinear  appearance  in  Fig.  6. 

An  uncertainty  estimate  for  the  eight  values  that  compose  the 
parameter  set  is  made  as  follows.  A  replicate  175-point  bottom- 
loss  data  set  is  created  by  randomly  selecting  with  replacement 
from  the  original  pool  of  bottom-loss  estimates  plus  associated 
measurement  variance;  this  is  called  a  Bootstrap  sample  [17]. 
This  sample  data  set  is  then  subjected  to  the  same  objective 
function  and  evaluated  over  the  constrained  parameter  search 
space,  yielding  a  new  parameter  set  called  a  Bootstrap  replica¬ 
tion  or  estimate.  This  process  is  repeated  1000  times  to  generate 
1000  replications  from  which  statistics  can  be  generated.  The 
uncertainties  reported  in  Fig.  6  represent  two  standard  devia¬ 
tions  computed  from  these  replications.  Note  that  constraining 
the  parameter  search  space  can  have  the  effect  of  artificially  re¬ 
ducing  estimates  of  uncertainty  and  our  approach  does  not,  as 
also  is  the  case  with  many  other  inversion  schemes,  quantita¬ 
tively  address  this  delicate  issue. 

Comparisons  between  the  measured  bottom  loss  and  the 
model  predictions  using  the  best-fit  geoacoustic  parameters 
are  shown  in  Fig.  7.  Over  the  2-20-kHz  frequency  range,  the 
measured  data  are  in  reasonable  agreement  with  the  model 
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(b)  Grazing  angle  (deg.) 


Grazing  angle  (deg.) 


Fig.  7.  Measured  bottom  loss  as  a  function  of  grazing  angle  and  frequency 
compared  to  model  for  —  20  log  | /2|  based  on  sediment  layer  overlaying 
homogenous  half-space  and  using  the  best- fit  geoacoustic  parameters. 

predictions  based  on  the  two-layered  sediment  reflection  co¬ 
efficient  model.  Importantly,  the  significant  oscillations  as  a 
function  of  grazing  angle  seen  in  the  data  at  4  and  8  kHz  are 
captured  by  the  model,  although  somewhat  less  so  at  8  kHz.  The 
oscillations  disappear  in  the  2-kHz  data  and  model,  presumably 
because  the  0.9-m  sediment  layer  becomes  nearly  transparent 
to  the  0.8-m  wavelength  at  2  kHz.  An  apparent  critical  angle 
given  by  arccos(ci/c3)  %  22°  is  suggested  by  both  the  2-kHz 
model  and  data,  with  this  angle  likely  applicable  to  frequencies 
<2  kHz.  Oscillations  also  disappear  in  both  16  and  20  kHz  data 
and  model;  in  this  case,  we  postulate  that  sediment  attenuation 
significantly  reduces  the  sound  intensity  reaching  the  interface 
between  the  sediment  layer  and  half-space  and  a  much  lower 
critical  angle  arccos  (ci/c2(0))  %  13°  appears  applicable  to 
higher  frequencies. 

Finally,  in  order  to  verify  the  assumptions  concerning  both 
plane-wave  reflection  and  ray-based  transmission  loss  calcula¬ 
tions  within  the  parameter  ranges  exercised  by  the  model,  we 
have  completely  duplicated  our  basic  technique  using  the  PE 
code,  based  on  the  range-dependent  acoustic  model  (RAM;  see 
[18]).  A  simulation  of  the  received  pulse  in  time  domain  is  ob¬ 
tained  via  a  Fourier  synthesis  of  the  narrow-band  complex  PE 
field  [19].  The  RAM  code  is  run  with  a  vertical-depth  discretiza¬ 
tion  of  our  bottom  model  shown  in  Fig.  6.  For  the  water  column. 


(b)  Grazing  Angle  (deg.) 

Fig.  8.  Bottom  loss  versus  grazing  angle  for  (a)  2  and  (b)  4  kHz.  The  solid  line 
is  the  predicted  value  of  plane  wave  reflection  model  and  circles  are  the  bottom 
losses  obtained  from  time-series  signals  simulated  using  a  parabolic  equation 
(PE)  code  with  the  geoacoutic  model  shown  in  Fig.  6. 

the  PE  is  subjected  to  vertical-depth  discretization  of  an  average 
of  several  CTD  casts.  The  procedure  involving  the  PE  is  per¬ 
formed  only  at  2  and  4  kHz,  as  assumptions  verified  at  these 
frequencies  will  hold  for  the  higher  frequencies.  Results  for  2 
and  4  kHz  are  shown  in  Fig.  8;  for  both,  the  predicted  values  of 
plane-wave  reflection  model  are  in  reasonable  agreement  with 
the  bottom  losses  obtained  from  the  time-domain  signal  calcu¬ 
lated  using  the  PE  method.  This  result  provides  the  validity  for 
both  the  ray  theory  and  the  plane-wave  assumption  over  the  pa¬ 
rameter  ranges  we  used  in  this  work. 

VI.  Summary  and  Conclusion 

Estimates  of  bottom  loss  made  within  the  frequency  range  of 
2-20  kHz,  and  within  the  grazing  angle  range  of  15°-24°  have 
been  presented  and  interpreted  with  a  model  for  the  modulus  of 
the  plane  wave  reflection  coefficient.  The  measurements  were 
made  in  the  East  China  Sea,  350  km  east  of  Shanghai  in  waters 
105  m  deep  (29°  39’  N,  126°  49’  E)  as  part  of  the  ASIAEX 
field  program  conducted  in  May-June  2001.  The  bottom-loss 
measurements  are  a  subset  of  propagation  measurements  made 
with  a  source  at  depth  25  or  50  m  and  two  receivers  (26  and 
52  m)  colocated  at  a  nominal  range  of  500  m.  Given  this  mea¬ 
surement  geometry,  along  with  a  bottom  reflection  rather  than 
a  bottom-scattering  modality,  these  measurements  represent  an 
inherent  horizontal  average  over  ranges  of  O  ( 10-100  m)  based 
on  the  dimension  of  the  first  Fresnel  zone. 

Geoacoustic  survey  measurements  [  1  ]  conducted  as  part  of 
the  ASIAEX  field  program  indicated  a  sediment  layer  of  O 
(1  m)  in  an  area  overlapping  the  experimental  site  and  coring 
data  showed  a  mean  grain  size  of  3.47  </>,  consistent  with  very 
fine  sand,  and  a  sound  speed  increasing  from  approximately 
1575  m/s  at  the  water-sediment  interface  to  1600-1675  m/s 
at  approximately  1  m  into  the  sediment.  The  core  data  in  [1] 
also  exhibit  considerable  vertical  variation  at  centimeter  scales 
or  fine-scale  layering.  In  situ  measurements  of  the  fine-scale 
bottom  relief  [2]  put  the  rms  bottom  relief  at  0.005  m  (within 
a  4-m  aperture). 

Knowledge  from  the  geoacoustic  surveys  motivated  the 
construction  of  a  forward  model  for  the  reflection  coefficient 
based  on  a  two-layered  sediment  for  which  the  sound  speed 
in  the  surficial  sediment  layer  is  allowed  to  vary  as  a  linear 
k 2  profile  and  the  region  below  this  layer  was  modeled  as 
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a  half-space  for  which  geoacoustic  parameters  are  fixed.  A 
set  of  eight  geoacoustic  parameters  was  identified  based  on 
a  weighted  least-squared  error  criterion  between  candidate 
parameter  sets  and  the  data,  in  an  exhaustive  search  of  the 
constrained  parameter  space  that  incorporated  a  priori  in¬ 
formation  about  the  seabed.  Estimates  of  uncertainty  for  the 
geoacoustic  parameters  were  calculated  using  the  Bootstrap 
method.  The  thickness,  density,  and  attenuation  of  surficial 
sediment  layer  were  estimated  as  0.9  ±  0.5  m,  2.0  ±  0.1  g/cm3, 
and  0.25  db  0.05  dB/m/kHz,  respectively.  Sediment  sound 
speed  in  the  surficial  layer  increased  from  1557  ±  4  m/s  at  the 
water-sediment  interface  to  1625  ±  35  m/s  at  depth  0.9  m.  The 
sound  speed,  density,  and  attenuation  coefficient  below  surficial 
layer  (or  in  the  half-space)  were  estimated  as  1635  ±  52  m/s, 
2.0  ±0.1  g/cm3,  and  0.25  ±  0.15  dB/m/kHz,  respectively. 
This  parameter  set  produced  model  curves  of  the  modulus  of 
the  plane  wave  reflection  coefficient  that  agreed  reasonably 
well  with  the  observations  of  bottom  loss.  The  parameter  set 
also  is  consistent  with  the  range  of  independently  measured 
geoacoustic  variables  discussed  in  [1].  Importantly,  the  model 
curves  reproduced  the  appearance  exhibited  by  the  data  of 
strong  variation  in  bottom  loss  as  function  of  grazing  angle  at 
4  and  8  kHz,  combined  with  reduction  of  such  variation  at  2, 
16,  and  20  kHz.  This  variation  in  loss  (reflection)  with  grazing 
angle  as  a  function  of  frequency  is  associated  with  0(1  m) 
-thick  layered  structure  in  the  near-surface  sediment,  which 
is  resolvable  by  our  mid-to-high-frequency  acoustic  measure¬ 
ments.  Our  inversion  results  apply  to  sediment  structure  at 
depths  less  than  about  3  m,  owing  to  in  part  the  combination 
of  lowest  frequency  being  2  kHz  and  maximum  grazing  angle 
being  approximately  25°.  Thus,  a  self-consistent  construction 
of  a  geoacoustic  model  for  this  area  should  necessarily  amalga¬ 
mate  the  mid-to-high-frequency  results  given  here  with  results 
obtained  at  lower  frequencies. 
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